CRITICAL VALUES OF RANDOM ANALYTIC FUNCTIONS ON 

COMPLEX MANIFOLDS 
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Abstract. We study the asymptotic distribution of critical values of random 
holomorphic sections s„ S H^(M™' , L") of powers of a positive line bundle 
{L,h) {M,u)) on a general Kahler manifold of dimension m. By critical 
value is meant the value of |s(2)|(,n at a critical point where S/fiS„{z) = 0, 
where Vj, is the Chern connection. The distribution of critical values of s„ is 
its empirical measure. Two main ensembles are considered: (i) the normalized 
Gaussian ensembles so that E||s„||^2 = 1 ^^'^ (ii) the spherical ensemble 
defined by Haar measure on the unit sphere SH^{M,L^) C H^(M, L") with 
||s„||^2 = 1- The main result is that the expected distributions of critical 
values in both the normalized Gaussian ensemble and the spherical ensemble 
tend to the same universal limit as n — > oo, given explicitly as an integral over 
m X m symmetric matrices. 



1. Introduction 

The purpose of this article is to determine the asymptotic distribution of critical 
values of random holomorphic polynomials of large degree, and their generalizations 
('holomorphic sections') to all compact Kahler manifolds (M, w) of any dimension 
m. We work in the same general setting as the articles [DSZl) [DSZ21 |DSZ3| on the 
distribution of critical points on Kahler manifolds and recall the definitions in ^ 
We first consider the critical value distribution of Gaussian random 'polynomials' s„ 
as the degree n — > cx) and then consider the more difficult and interesting problem 
of critical values of normalized random polynomials with HsuHl^ = 1- We refer 
to the latter as the spherical critical value distribution since the 'polynomials' are 
drawn at random from the unit sphere in Hilbert space. We regard the spherical 
distribution as primary for the critical value distribution since a critical value is 
only counted once in a line {cs„,c G C} of sections and one can relate the heights 
of the critical values to the other threshold heights of normalized 'polynomials'. 
Theorem [l] shows that a special normalized Gaussian critical value distribution has 
a universal limit independent of the manifold and Kahler metric. The main result 
of this article, Theorem [2j shows that the spherical critical value distribution has 
the same universal limit. We also give a limit formula for the simpler spherical 
value distribution in Theorem [3j The spherical limit results may be viewed as a 
Poincare-Borel theorems for critical values. 

The spherical distribution of critical values is potentially useful in analyzing 
the Morse theory of the modulus |s„(z)|^„ of random sections, since it is at critical 
values that the level sets change topology. They are also important in understanding 
the topography of random surfaces, i.e. the graphs of the modulus y — |s„(z)|^„ 
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of 'polynomials' of degree n, which are often visualized as mountain landscapes 
above some given sea level. It is known that sup^.^^^,/ \\sn{z)\\h" < Cm'^l'^ when 
IIskIIl^ = 1 and dime Af = m and it is proved in [SZ^ (after a long history of 
similar results in other settings) that the expected sup norm of such normalized 
sections is bounded by a universal constant times \/\o%n. Thus in a measure sense, 
typical 'polynomials' of degree n and norm one have global maxima < C-\/logn, and 
conjecturally the median should be of this form for some C Deterministically, 
critical values of all normalized polynomials lie in [0, Cn™l'^\ . It would be interesting 
to know the exact height (for I? normalized polynomials) at which the peak of 
the random mountain first occurs. At a certain threshold height, the mountain 
tops are sometimes conjectured to form a Poisson spatial process, and it would be 
interesting to know the connectivity properties of the landscape at lower sea levels. 
The calculation of the spherical density of critical points is only a calculation but 
it is probably a necessary one for the more involved landacape questions. The 
Kahler setting is a model for other settings in which one studies normed random 
waves with a notion of degree or eigenvalue, such as random spherical harmonics or 
more general Riemannian waves on Riemannian manifolds, in which the principal 
modification is in the asymptotics of the relevant covariance functions. 

Before stating the results, we introduce some notation and background. Com- 
pact complex manifolds have no non-constant holomorphic functions and the nat- 
ural replacement for them are twisted holomorphic functions know as holomorphic 
sections of complex line bundles tt : L — >■ M. Here, the fiber over z € M is 
a one-complex dimensional space and a holomorphic section is a map s : M ^ L 
satisfying ds = 0,7r o s{z) = z. Degree n sections are sections of the nth tensor 
power L" of L, and the space of holomorphic sections is denoted 7J"(M, L"). Its 
dimension is given asymptotically by 



When dime M — 1, i.e. when M is a Riemann surface, then the natural examples 
are polynomials of degree n {g = 0), theta functions of degree b {g — I) and 
holomorphic differentials of type (dz)" for g > 2. The techniques and results of this 
article, as in the predecessors |BSZ1 IDSZI) . hold in this general geometric setting. 

The Kahler metric u determines a Hermitian metric h and connection V on i 
and on its powers. The Hermitian metric satisfies 99 log h = oj, and the connection 
V is known as the Chern connection and is compatible with the Hermitian metric 
h on and complex structure on L. As recalled in f|2j the Hermitian metric h and 
Kahler form lo give rise to a definition of Gaussian random holmorphic section in 
i?0(M, L"). 

By a critical point of a holomorphic section s e H^{M,L'"-) we mean a point 
z e M so that 



Thus the section is 'parallel' at z. Equivalently, critical points are points where the 
norm square is critical d|s„|^„ = and so we are studying the critical points and 
values of the real-valued function |s(z)|^j„. More precisely. 



(1) 



dime H 



(A/, L") ~ Cn 



(2) 



Vs(z) = 0. 



(3) 




The methods of this article give an explicit value of C, which we defer to a future article. 
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We note that a; = is surely a critical value of |s„(z)||„ since every section has 
zeros. However we omit the zero value in the definition of the empirical measure 
of critical values. When the Hermitian line bundle (L, h) is positive, as we assume, 



the only local minima of \s\h are its zeros. Therefore the critical points in (13) are 
either saddle points or local maxima. 

To define Gaussian random holomorphic sections, introduce a family of Gaussian 



measures adapted to the Hermitian metric and the associated inner product ( 20 ) 
on sections. For any a > we put 

(4) d7a(s«) = i^r-e'^\^\'da , s„ = £ a^s] , 

where {s",--- , sJJ } is the orthonormal basis of H'^{M,L'^) with respect to the 
inner product (20). Equivalently, the coefficients a" arc complex Gaussian random 



variables which satisfy the following normalization conditions, 

(5) E>^ = 0, E>^7 = ^4,, E»;' = 

Here, we denote the expectation with respect to 7^ by E ^ . Under this normaliza- 
tion, we have the expected norm of s„, 

(6) E,^||.„||^„ = ^. 

a 

Thus the covariance kernel of 7^ is 

(7) ^l{s{z)J{w)) = -n,,{z,w), 

a 



where H„ is the Szego projector with respect to (20). 

When a = d„ we call the ensemble the normalized Gaussian measure. It is given 

by 

(8) dp^l{s^) = (^)'i"e-'^'>l'da, 

TT 

and from (|6| it follows that 

(9) EM|s„||2„ = 1. 

As will be seen below, the density of critical values has a limit for this sequence of 
probability measures. We discuss the relations of the densities as a varies in f|3] 

The distribution of critical points of a section is defined by the un-normalized 
empirical measure 

(10) - 

z:Vs(z)=0 

The Chern connection V associated with h is not holomorphic, and the number of 
critical points depends on the section s. The statistics of the number of critical 
points in the normalized Gaussian ensemble was determined in |DSZ1[ IDSZ2| . The 
critical point distribution is invariant under s — )■ cs and therefore it is equivalent 
to work with Gaussian or spherical distributions. It is proved in Corollary 5 of 
|DSZ2| that the expected number J^^qh critical points of Morse index q for any 
positively curved Hermitian metric ft, on a Kahler manifold of dimension m satisfies 



(11) AC™* 
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Here, 6o<j is a Betti number and ci{L) is the first Chern class, both of which are 
topological invariants. For instance, in dimension one there are roughly |n saddles 
points and |n local maximal for Gaussian random polynomials of degree n with 
the SU (2) (or Fubini-Study) inner product. 

1.1. Statement of results. In this article we study the Gaussian, resp. spherical, 
distribution of critical values, 

(12) CVs := {\siz)\h^ : Vs{z) = 0, s G i/"(M, L")} 

in the limit as the degree n tends to oo. Thus, the "value" of the section at a critical 
point is the Hermitian norm in K of the section; since s{z) g it would not make 



sense to study the values themselves. In view of (11), we define the (normalized) 



empirical measure of nonvanishing critical values of |s„|^n by 
(13) ^^^ = ^^ E ^l-l-- 

z: V's„{z)=0 

Note that it is not necessarily a probability measure but from the results of |DSZ2] 



(such as (11)) it follows that for any e > there exists a constant C so that 
#{z : Vs(z) = 0} < Cn™ except for a set of sections of measure < e. 

We define the Gaussian density of critical values D"(a;) as the expected density 
of E " CVs in the sense of distribution, 

(14) E^(4^ J2 ^i\'n\h-)]^[ i^{xmx)dx for i, e CoiR+) 



where dx is the Lebesgue measure on M. In ^we will calculate the densities for all 
a. The Kac-Rice are quite complicated for fixed n, but for the normalized Gaussian 
ensemble there are simple asymptotics. 

To state the result we need some notations. We denote by 5'(C'") = C 2 the 
space of complex symmetric matrices. We also denote by the Legesgue measure 
on S'(C™). We also define the special matrix P {Sjj'Sqq' + Sjq'Sqj') ^ 



Theorem 1. Let {M,uj) be an m- dimensional compact Kdhler manifold. Let 
{L, h) — >■ M be a polarized positive holomorphic line bundle. Let DJ^" be the expected 



density of critical values defined in {I4) with a = dn, i-e. the normalized Gaussian 
density. Then we have, 

(15) Bi"{x) = B^{x) + o(J^{x{l + x^)e~^")^ on (0,oo), 

where 

T>oc{x) - /,„(|x|)|x|e-l^l', with frn{\x\) = c„, / e'l^l' IVpcP - d^- 

JS(C'") 

Here, c„i ~ (m+iKm+2) where V is the volume of {M,uj). The asymptotics can 

TT 2 

be differentiated any number of times ( with appropriate changes in the polynomial 
growth in the remainder estimate.) 



Remark 1. Henceforth we generally set V = 1 for notational simplicity. 
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In the case of Riemann surfaces when to = 1, we have P = 2. Assuming the 
volume of M is tt, then 

(16) D^"(x) =x(2x2-4 + 8e--j e^" +0(^ - 

Below is the computer graphic of the leading term, 




Graph of Doo{x) in dimension one ^ 

The critical point densities for D" have the scaling relation D"(a;) = ai'Dl^{a^ x) 
(Lemma[2]) , and from this one can determine the asymptotics for the other Gaussian 
ensembles. We also give a similar formula for the simpler expected distribution 
function of random sections in f|9] 

The proof of Theorem [T] is is based on the Kac-Rice formula in Lemmas [3] and |4j 
which give exact formulae for all of the D"(a;). The asymptotics then follow from 
the complete asymptotic expansion of the covariance kernel Q in S 7.1 In the case 
of SU{m +1) polynomials on CP™ we give an exact formula for all rt in f|6] 

Remark 2. As in [DSZ2j Theorem 1.2 (see also |Bauj ). we can give similar for- 
mulae for the distribution of critical values when the critical point is constrained to 
have a specified Morse index. The formula only changes in that we integrate over 
the subset S'^{C"^) of matrices of index q. 

1.2. Density of critical values in the spherical ensemble. As mentioned 
above, the critical point distribution is homogeneous, i.e the same for s and cs. 
The critical value distribution is however not homogeneous, since the critical values 
are multiplied by |c|. Since the mass of the normalized Gaussian measure is asymp- 
totically concentrated near the unit sphere as d„ ^ oo, the critical values of a line 
of sections {cs} are weighted most for values of |c| close to 1. This weighting is a 
re-scaled version of the one in the classical Poincare Borel theorem, which states 
that the spherical probability measure i^d on the sphere S'^{y/d) tends to the Gauss- 
ian measure as d — oo. More precisely, if Pd ^ M*^' is Pd{x) = \/d{xi, . . . , Xk), 
then for aU k, Pd*Vd ^ Ik = {2'K)-'^/'^e-\^\^ /"^dx . Moreover, 

7d{x e W : ||a;||2 > < e-'"''/\ 7d{a; G M'' : \\x\\^ < (1 - e)d} < e"^'^/*. 
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Our spherical probability measure is normalized Haar measure dvn on 

(17) SH°{M,L'') = {se H''{M,L'^) : \\s\\lo_ ^ 1}. 

We refer to the corresponding probability space as the spherical ensemble. What 
we are calling the normalized Gaussian measure ([s]) concentrates exponentially on 
this unit sphere. We denote the expectation with respect to dvn by E and define 
the (normalized) spherical density of critical points D^(a;) by 

(18) Y^y^CVs^'Dn{x)dx. 

In fact it makes more sense to pass to the quotient Fubini-Study probability measure 
on the projective space PiJ*'(Af, L^) of sections since the critical value distribution 
is invariant under multiplication by e*^. 

We view the spherical density ([Ts]) as primary because fixing ||s„||i2 = 1 sets 
a scale against which one can calibrate the heights at which interesting features 
of the landscapes occur. For instance, as mentioned above, the spherical critical 
value distribution is supported in [0, C\/c?„], and its median should occur at a 
constant times -^/log n. In the sequel we plan to study such distinguished levels in 
more detail. The main result of this article is: 

Theorem 2. The density of critical values in the spherical ensemble SH^{M,L") 
on any compact Kdhler manifold has the universal limit, 

lim T)l{u) = D^{u). 

n— >-C30 

Thus, the spherical critical point distribution tends to the same universal limit as 
the normalized Gaussian measure of Theorem [T] As the asymptotics (and Graph) 
indicate, the most probable critical value and the median of the critical value dis- 
tribution is around 1 in dimension one. An upper bound for the median may 
be derived from the exact formula for the spherical critical point density. In a 
subsequent article we will apply Theorem [2] or more precisely its proof to obtain a 
formula for the median and for the asymptotics of the critical point density with 
fixed Morse index in special n-dependent intervals. 

As discussed above. Theorem [2] is a Poincare-Borel type theorem for critical 
values. Intuitively it is based on the concentration of normalized Gaussian measure 
around the unit sphere, but in its details it uses the special scalings of the critical 
value distribution and the asymptotics of the covariance kernels and therefore does 
not seem to follow directly from the classical Poincare-Borel theorem. The proof 
is based on a Laplace transform relation between the spherical and normalized 
Gaussian critical value distributions. 

The Poincare-Borel relation between the normalized Gaussian and spherical ex- 
pectations of the critical value distribution holds also for the full value distribution. 
We denote the density of values in the spherical ensemble by f^{u)du. In ^we 
prove: 

Theorem 3. The density of values in the spherical ensemble SH'^{M,L") of any 
Kdhler manifold has the limit, 



lim fn{u) = 2ue 



CRITICAL VALUES 



7 



1.3. Related results and problems. Other articles over the last ten years de- 
voted to the statistics of critical points of Gaussian random fields include |DSZ11 
IDSZ21 IDSZ31 iBiHl iMiHl iNSVl E [GWl EBXI IHl [Z] . The Kac-Rice formula 
for the Gaussian critical value distribution of holomorphic sections were originally 
obtained in |SZ3) but the asymptotics were not determined as explicitly as in this 
article. As mentioned above, we view the Gaussian formalism mainly as a method 
for computing the spherical distribution. In the real domain, a Kac-Rice formula 
and an asymptotic analysis of the critical point distribution are given in jN2| . 



The expected value of the empirical measure ( 13 ) is only the first and simplest of 
the many probabilistic problems on critical values. As in the case of zeros or critical 
points, one may ask for the variance of linear statistics (pairings of smooth test 
functions with the empirical measure), the asymptotic normality of linear statistics, 
large deviations properties and so on. As mentioned above, the precise structure 
of the landscapes defined by ?; = |sn(2;)|/i" is unknown in many respects. 

In the case of polynomials of one complex variable p{z), one might instead use 
the standard complex derivative ^ to define critical points and critical values, but 
it is in fact a meromorphic connection with a pole at infinity and leads to quite a 
different theory. To our knowledge, statistics of critical points in the latter sense 
have only been studied in [HI IFW| . In [FW, , the authors studied the expected 
density of critical values of Gaussian SU (2) random polynomials p„ defined on C 
with respect to the meromorphic connection Also, in TT, B. Hanin studies the 
correlation between zeros and critical points for this classsical connection. 

The main result of |FWj is that the (un-normalized) expected density of nonva- 
nishing critical values of the modulus of \pn \ satisfies 

V^: |P.|'=0 / 

on as x away from where \pn\' = ^g^"^ . This result is quite different from 
Theorem [l] due to the fact that the connection d/dz is a flat meromorphic one 
rather than the smooth but non-holomorphic connection of this article. 



2. Background on Gaussian measures in the Kahler setting 

2.1. Kahler manifolds. Let {M,ijj) be a compact Kahler manifold of complex 
dimension m with the local Kahler potential w = ddtp. Let (L, h) ^ M he a 
positive holomorphic line bundle such that the curvature of the Hermitian metric 
h, 

(19) Qh = -dd\og\e\l 

is a positive (1,1) form GHJ. Here, e is a local non- vanishing holomorphic section 
of L over an open set U C M such that locally L\u ^ U x C and |e|/j = h{e, e)^/^ 
is the pointwise /i-norm of e. 

We denote by H^{M,U'^) the space of global holomorphic sections of L" — 
L ® ■ ■ ■ ® L. Under the local coordinate, we can write the global holomorphic 
section as s„ = fn^®^ where /„ is a holomorphic function on U . We denote the 
dimension of iJ° (M, L") by d„. 
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The Hermitian metric h induces a Hermitian metric /i" on i" given by le**"!/!" — 
|e|5J. Throughout the article, we assume the polarized condition 0^ = uj such that in 
the local coordinate, we have the /i-norm \e\h = and hence |s„|ftn = |/n|e~~f . 

We decompose the Chern connection V = V' + V" of the Hermitian line bundle 
(L", /i") into holomorphic and antiholomorphic parts where in the local coordinate 
V' = 4 + ndlog h and V" = dg [GH] , 

We can define an inner product on H'^{M, L") as, 

(20) = / h-{s^,sl^)dV 

Jm 

where dV = ^ is the volume form. We recall that throughout the article we 
assume the volume is normalized as dV = 1. 
We write this in the local coordinates, 

(21) - / hhe-^^'dV 

JM 

where = /fe®" and s^' = /a^e®". 



2.2. Gaussian measures. Recall that a Gaussian measure on R" is a measure of 
the form 

7A — ; — , dxi • • • dxn , 

' (27r)"/2Vd^ 

where A is a positive definite symmetric n x n matrix. The matrix A gives the 
second moments of 7a: 

(22) {xjXk)j^ = Ajfc . 

This Gaussian measure is also characterized by its Fourier transform 

(23) 7^(ti,...,t„) =e-5i:A,.*.*. . 

If we let A be the n x n identity matrix, we obtain the standard Gaussian measure 
on M", 

with the property that the Xj are independent Gaussian variables with mean and 
variance 1. 

A complex Gaussian measure on C'^ is a measure of the form 

^- = ^Me^^'^' 

where dz denotes Lebesgue measure on C'^, and A is a positive definite Hermitian 
k X k matrix. The matrix A — (Aq,^) is the covariance matrix of 7a : 

(25) ( = A„^, l<a,/3<fc. 
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3. A ONE-PARAMETER FAMILY OF GAUSSIAN MEASURES AND THEIR CRITICAL 

POINT DISTRIBUTIONS 

The one-parameter family of complex Gaussian measures Q on i/°(M, L") may 
be written formally as 

TT 

where Ds is Lebesgue measure. 

The normalization comes from the calculation 

(26) 

where ujk = r{k/2) ^^'^ surface measure of the unit sphere S'^^^ C M'^'. 

As in ( 14 1 we denote the normalized density of critical values with respect to 7" 

by 

(27) E^Cy, = DJJ(x)dx. 

We note that D"dx is not a probability measure on R+ since CVs is not in 
general a probability measure. The mass of the measure can be determined from 
the eventual Kac-Rice formulae of SjS] 

3.1. Normalized Hemitian Gaussian measure. As mentioned in the introduc- 
tion, when a = dn it's the normalized Gaussian measure ([s]) which is characterized 

by 

(28) Ea^ = 0, Ea^5' = J-4,, Ea^J^O 
Under this normalization, the expected norm of s„ is 1 ([T]) . 

3.2. Relations between Gaussian critical value densities. Next we compare 
densities (27 1 as a changes. The first step is 

Lemma 1. For any s G H'^{M,L"') with non- degenerate critical points, and for 
any r > 0, 

CVrs = CVsir-^x). 
Proof. For any / e C(M) we have, 

{CVrs,f)= J2 fir\siz)\h^)^{CVsix)J{rx)). 

z:Vs(z)=0 

Changing variables to y = ra; completes the proof. 

□ 

Since D'^dx transforms by the inverse dilation, the density has the transforma- 
tion law, 



Lemma 2. D"(a;) = D^(q2 x). 
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Proof. We have, 

where we change variables s — > s' = q^^s and apply Lemma [l] □ 

Of course this implies that 'D"{x)dx all have the same mass. Since a — d„ in 
the normalized Gaussian measure, we have: 

Corollary 1. T)i-{x) = dk r>i{dkx) 

4. Kac-Rice formula for the Gaussian critical point density 

The main result of this section (Lemma |3]) gives a genric Kac-Rice formula for 
the expected density DJ^" of critical points with respect to the normalized Gaussian 
measure. Very general Kac-Rice formulae applicable to the critical point density in 
our setting are given in |BSZ| . Other presentations can be found in |ATi IN1| . The 
historical references are [Kl |R] . 

The Kac-Rice formula is as follows, let f{t) be a real valued smooth stochastic 
process on the finite interval / C M. Then the expected number zeros. 



(29) 



E #{<€/: f{t) = 0}= f I \y\pt{0.y)dxdt 



where pt{Q,y) is the joint density of {f,f') evaluated at (0,y), dy and dt are 
Lebesgue measures on M. If / is a Gaussian process, then the joint density pt(x, y) 
is uniquely determined by the covariance matrix of (/, /') [BSZl ET] . 

4.1. Kac-Rice formula for the critical point density. In this subsection, we 
will derive the Kac-Rice formula for the expected density DJ^" of critical values of 
jsnU" with respct to 7^". As we will show this particular Gaussian density has 
a limit as n — > (X). The formula may be derived from |BSZ1 [DSZl] but we take 
advantage of some simplifications to speed up the proof. To simplify notation, we 
write 

D — D"*" 

■— '-'n ■ 

In the local coordinate U = C™ and a local trivialization of L, we write the 
normalized Gaussian random sections as, 

(30) s„ = /„e«". 
where 

(31) /» = Eaj/; 

where {uj} are normalized Gaussian random variable ( 28 ) and locally {s" = /"e**", j = 
1, • • • ,dn} is an orthogonal basis of H^{M,L'^) with respect to the inner product 

(po|. 
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The smooth Chern connection then has the form [GH| 
(32) V's„ = (/; - n9^/„)e®". 

Here, d = d + d is the decomposition into terms of type (1, 0), resp. (0, 1). Here, 

- 9/ 
■ij dzj 



fn — 9f = J2j ^dzj. Thus, in the local coordinate, the empirical measure (13) 



has the form 

(33) CV, = ^.Y.^ 



where 

(34) r! = {zeC": (/;-na^/„)e-^ =0}. 

We also introduce the locally defined empirical measure of complex critical values 

(35) CVs = — y<5, 

We will determine the expected density D„ of CY s on C and then integrate out the 
angle variable D„(|a;|, 0)|a;|d0 to obtain the (global) expected density of GVs- 
In other words, we use that (V', CV^ — (V', CV ^ for radial functions and thus, 
CVs — -K^CV s where tt : C* M+ the map z — > \z\. 
The result is: 

Lemma 3. The expected distribution of critical values CVg is given by the formula, 

(36) D„(x) - 4^ / / / ^' 0' |det - n^Ix^) \ xd^dYdO 

n Jo J M J s(C"^) 

where ^ € ^(C™) = C ' 2^ ' is the space ofmxm complex symmetric matrices and 
d£, is the Lebesgue measure on S{C"^), p"(a;,0, 0,^) is the joint density o/p"(t/,0,^) 
of normalized Gaussian random variables {fm fni fn) evaluate at f'^ — 0, here we 
substitute y := xe*^ by the map n. (The formula of p"{y,0, is given explicitly in 
Lemma^. 

Proof. We first introduce some notations: 

Pn = fne~^ e C, g„ = (/; - nd^f„)e~'^ G C™, r„ = /^'e""^ G S{C), 

then pn, and r„ are all complex Gaussian random variables. 

By definition of the delta function, we have for any test functions ip G C^(C), 



n™ ^ 2 ' ^' 



n" 

z: q„(z)=0 



1 

1 



So{qn)^iPn)dqn A 



^o(gn)'0(p«) 



det 



%i|2 _ 1^12 
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where dz is the Lebesgue measure on C. By direct computations, we have, 



and 



— = -nddippn - -dipq^ 
oz 2 



where and are mx m symmetric matrices; for simphcity, for any matrix A 
we denote, 

|Ap -.^ AA* 
By taking expectation on both sides, we have. 



2eo 



n" 
1 



det 



dqn ,2 _ |^gn|2 



dz 



i^{y)Pz{y. 0, |det (le + n^^^a,^*^ - n^VyP - n'\ddip\''\y\'') \ d^dzdy 



where Pz(y, s,^) is the joint probabihty of the Gaussian random field {pmqmfn) 
and dx is the Legesgue measure on C. 

Thus, the expected density E CV s is, 
(37) 

D„(2/) = 4^/ / p,Xy,Q,0\'^^^{\i + n'dvd^*y-nd\y\''-n'\dd^\''\y\^)\d^dz 
We rewrite y in polar coordinate {x, 9). Then the expected density of CVg is, 

(38) D„(a;) = / f^n{x,e)xd9 

Ja 

The next step is to get an explicit geometric formula for D„(x). The empirical 
measure CVs and its average D„ are independent of coordinates and frames on the 
Kahler manifold M and line bundle L". We may choose Kahler normal coordinate 
to simplify the above integral. 

We freeze at a point zq as the origin of the coordinate patch to simplify the 
integrand at zq. It is well known that in terms of Kahler normal coordinates {zj}, 
the Kahler potential ip has the expansion in the neighborhood of zq: 

(39) (p{z, z) = \\zf Rfkpq{za)zjZ-kZpZq + 0(||2;||^) . 

In general, (f contains a pluriharmonic term f{z) + /(z), but a change of frame for 
L eliminates that term up to fourth order. Thus 

(40) dLp{zQ) = 0, aV(^o) = 0, ddip{zo) = 1, dV{zo) = dz. 

Such frames are called adapted in |BSZ[|DS^ . Hence, the joint density of r„) 
at Zq is the same as the joint density of Gaussian process (/„, /^'). 
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Thus by (|37])(38)(40), we obtain the global expression, 



r>n{x) = ^f f [ p:{x,e,0,0\det{\^f~n^Ix^)\xd^dVd9 
n Jo Jm Js(C'") 

which completes the proof. □ 

5. Calculation of the joint probability density 

In this section we calculation the joint probability distribution p" of Lemma [3] 
with respect to the normalized Gaussian measure 7^" . 

5.1. Density (a;, s, ^). In this subsection, we will derive the formula for s, ^) 
of the joint density of the Gaussian process (/«, /4, /")• It is given by the formula 

(41) p:iy^s,o = id^e-p<! -( I - 1 '(A")-^ I 

where jg the dimension of the Gaussian process {fn, fn: fn) s-'^d A" is 

the covariance matrix of this process. 

We rearrange the order the Gaussian process and write A" as the covariance 
matrix of (/^, /„), then we rewrite, 





(42) P?(.,.,0 = i^exp 




The covariance kernel is defined by 

(43) n„(z,«;):=£/;(z)/,"H. 

where {/"} is defined in (30 1. 
Then, we have 

(44) E(/„(z)/„(u;)) = -^n„(z,ii;) 

The notation Iln{z,w) usually refers to the Szego kernel but in fact (44 1 is the 
Bergman kernel, which has the pointwise TYZ expansion [Ll ITl [Ze] . 

(45) n„(z, z) = n™e"^(^) [l + ai(z)n-i + 02(2)71-^ + • • • ] , 

where ai is the scalar curvature. Integrating over M with respect to e~"''^dV gives 
the well-known dimension polynomial, 

(46) d„ = n"(l + n"M aidF + n^M a2dV + ---) 

Jm Jm 

The covariance matrix is then given by. 



where 

(48) 

dn ozaw 
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is an TO X m matrix; 



(49) 



IS an TO X 



Br, — 



matrix; 



dzd'^w \^=w 



(50) 



Cn — 



1 



^ m +m+2 ^ m +m+2 



matrix. 



d^n„(z,w) I 

g2^ \z=w 



dz U='^ 



^n„{z,w) I 



n„(z,2:) 



Thus, we have 
Lemma 4. Wif/i f/ie above notations, 

1 1 



(51) 

where 
(52) 



tt'''" det A„ det A„ 



exp 



A, 



"(Cn — B*^Aj^^Bn). 



6. Calculation in the Fubini-Study case 



In this section, we give exphcit formulae for the Kac-Rice density of critical 
points for SU{m + 1) polynomials. This is the case where M = CP™, where L is 
the line bundle 0(1) whose sections are linear functions on C™+^, and so sections 
L" = 0{n) are homogeneous polynomials of degree n. We equip 0(1) with its 
Fubini-Study metric /ifs given by 

Uw)\ 



(53) 



|s|i/lFs(N) = 



\w\ 



w = {wq, . . .,Wm) e C 



m+1 



for s e C'"+i* = H"{CV"\0{1)), where \w\^ = \wj\^ and [w] e CP™ denotes 

the complex line through w. The Kahler form on CP™ is the Fubini-Study form 



(54) 



'^Fs = ©/iFs = ^^<991og|w|' 



We denote by D 
Proposition 1. 

(55) Df^(™+i)(.T)=c,„a;e 



SU{m+i) density of critical points for 7^". Our result is 



5(C'") 



det 



n- 1, 



ere c,„ is as m Theorem 



SU{m+l) 



We note that the only difference between D 
rem [1] is in the limit — >■ 1. In dimesnion to = 1, 



and the limit Doo of Theo- 



(56) 



nsu{2)( 



X ) — —xe 



n-1 



'-2r- Ix' 



dr. 
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Proof. The Szego kernel for 0{n) is given in an affine chart Zq = 1 with Zj = ^ 

by 

77-1-1 

nffO(cpi,o(„))(^,w;) = ^(1 + zw^e^^iz) e®"(«;) . 

Since our formula is invariant when the Szego kernel is multiplied by a constant, 
we can replace the above by the normalized Szego kernel 



(57) Ur,{z,w):={l + zw)'' 

^{z) = -\og\e{z)\ls=iog{l + \z\^) 



in our computation. 
Since 



we have 



^(0) = ^(0) = g(0) = 0; 



dz dz"^ 

i.e., e„ is an adapted frame at z = 0. Hence when computing the (normalized) 
matrices -B„, C„ for i?°(CP^, 0(n)), we can take the usual derivatives of Iljv. 
Indeed, we have 



dz 
d'^z 

dzdw 
dzd'^w 
d'^zd'^w 



It follows that 



n(l + zw)^ , 

n(n- l)(l + ^;u))"-^(u))2 

n{n-l)il + zwr-^z'^ 

n{l + zw)"~'^ + n{n - 1)(1 + zw)'^~'^zw . 

2n{n - 1)(1 + zw)"-^z + n{n - l)(n - 2)(1 + zw)"-^z^w 

2n{n - 1)(1 + 0w)"-2 + 4n(n - l)(n - 2)(1 + zw)''-^ziB 
+n(n - l)(n - 2)(n - 3)(1 + zw)''-'^z'^w'^ 



(58) A„ = = (nd.,)..^^ 

(59) ^" = ^(^^U=. ^|_)=(0 0) 
(an m x matrix); 

1 / d^n^jz.w) I a^n„(z,^;) , \ /2n(n-l)/ 0^ 

(60) C„ = ^ > ''^^ = 

dn['-^^U=^ U^{z,z) ) [ 1, 

(an '"^+^'"+^ X "'+^"+^ matrix.) 
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Since 5 = 0, 



dn 



n{n - l)P 
1 



where 
(61) 
Thus 
(62) 
and 
(63) 



det A„ = (n^)^^ det P, 



'(n(n^l)r^p-^\i\^-\y\^ 



Write y in the polar coordinate {x,6) and combine Lemmas [s] and [4] we have 



(64) 



|det(|C| 



1 xe 



2n 



2 2 2 

n X 



M J5(C™) 

) I didVdO 



1 

detp' 



Now we change variables ^ — >■ \Jn{n — 1)^, apply the assumption J^^ dT^ = 1 
and integrate Q variable to get, 



(65) Df^("+i)(a;) = 



27rxe ^ 



S(C-) 



detP 



det 



n- 1 



Now change variables ^ — > v^'C to obtain the stated result. 



□ 



For m = 1 we get 

Dr(^)(x) 



-xe Jo e 



^ ^ — r I o n— 1 



r — X \ dr 



(66) 



Jo 



2r 



n— 1 



2x2-4;^+8exp{-^t-i:r2} 



7. Proof of Theorem [T] 

In this section, we will complete the proof of Theorem [T] The main additional 
ingredient is the asymptotic expansion of the Szcgo kernel. Otherwise the compu- 
tation is similar to the case of CP™. 

7.1. Covariance matrix. We now calculate the leading terms in An and A„. The 
key point is to calculate the mixed derivatives of n„ on the diagonal. It is convenient 
to do the calculation in Kahler normal coordinates about a point zq in M. 



Now take two derivatives on both sides of ( 45 1 , 



(67) 



+ 23?aj(/?e"'^[9jvai + •••]+ n^dj(pdj>ipe'"^[l + n^^ai + • • • ]) 
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We apply (39) at the origin zq to get, 

Ai{zo) = djdfUnizo, zo) = — {nl + aj + n~^{a2l + d.jdj>ai + •••)) 
By the same arguments, we compute B{zq) up to the leading order term, 



(69) B„(zo) = 9j9/9g'n„(zo,2:o) = ^fe'C?g'ai+<5j9'5/aiH , n ^5jai + - 

and the leading terms in each entry of C„(zo), 
(70) 



C„(zo) 



n'^{5jj:5qq> A- 5jqi5qj'){l + n ^ai + ---) n ^djdqai+n ^djdqa2 + 



1 + n~^ai + 



n ^djdqai+n ^djdqa2 + 

We refer to jDSZ2j for more details of those computations. 

Together with Lemmas [3] and |4j the formulae for An , _B„ , C„ give an explicit 
formula for D„ on a general Kahler manifold. 

7.2. Proof of Theorem [ij In this section we calculate the leading order term as 
n — >■ oo of the expected density D„(a;). We have, 



dn 



(C„ — -B,ij4„ ^-B*) 



n^P + 0{n) n-^djdqai+0{n-^) 

dn {n-^djdqai + Oin-^) 1 + n-^ai + Oin-^) 

n^P + 0{n) n-^djdqai + 0{n-^) 

n-'^djdqai+0{n-^) 1 + n-^ai + 0(n-2) 



in the last step, we apply the full expansion of dn (46 ) and where P is given in (61 1. 
Thus 



(71) 



det 



(n^)"^ det P. 



Applying ( 46 1 to ( 48 1 again, we obtain the asymptotic analogoue of ( 62 ) , 
(72) 



det An ^ n" 



We then have, 
(73) 
Thus, 

(74) p.(y,0,O = 
hence by Lemma [3j 

T>nix) 



Ti P ^ n ^d d 

?^ f ^ + lower order 



1 



terms 



1 



detP n™n'"^+™ 



1 xe ^ 



(75) 



e " ^ + lower order terms. 



1 



detP 



'M ^S'(C'") 

|det (I^P - n^x^) I d^dVdO + lower order terms 



here we substitute y = xe in the expression of ( 74 ) . 
The remainder estimate is discussed in 37.31 below. 
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Now we change variables ^ — ^ n^, apply the assumption Jj^^ dV = 1 and integrate 
in the 6 variable to get, 

(76) D„(x) = / ;^e-^-^l«l^ |det (1^1^ - Ix') \ + 0{'-) 

Now change variable ^ — > y/P£_ to get, 

B„{x) = c„,xe--' /src-O det (j^/P^p - /x^) rf^ + 0(i) 

(77) 

= /™(x) + 0(i). 

7.3. Remainder estimate. To prove the remainder estimate we use more terms 
in the Bergman kernel expansion and Taylor expansions of the various functions in 
the Kac-Rice integrand. We sketch the proof as follows. First we have, 

A„ = — (C„ - BnA~ B*) 

_ /n2p(l_+ n-ifli + ...) n-^d^dqai + 0{n-^) 
~ d„ \n-^djdgai + 0{n-^) I + n-^a^ + 0{n-^) 

_ 1 /n2p(l + n-ifli + • • • ) n-^djdgai + 0{n-^) 

^ l + n-iai + --- \ n-^djdqai+0{n-'^) 1 + n-^ai + 0(n-2) 

Hence there exists a complete asymptotic expansion, 

det A„ = (n^)^— detP(l + +rr'^b2 H ) 

Same apply to 

det An = ^"(1 + n"^ci + 7i"^C2 H ) 

where bj and Cj are polynomials of curvature and uniformly bounded. 
This two identities imply the full expansion, 

det A = det A„ det A„ = n™'+^"(l + n^^di + n^^da + • • • ) det P 

We also have, 

1 _ /n-2p-i(l-n-iai + ...) ^-^(a.a.ai + • • • )A 
^'^^ " ~ 1^ n-^djdgai + ---) l-n-iai + --- . j 

It follows that 

Pziy^^iO = 7r<i>"(l+„-idi+n-irf2+--- )dctP „m„m2+™ 



where (, ) is the nature inner product on C. 

We substitute this formula into the following integration in Lemma [3] 

^n{x)^^^f f f p"Ax,9,0,O\det{^C -n'lx')\xd^dVd9, 
then we change variables ^ — >■ and y — xe^^ to obtain the exact formula, 

^ f f^"" f (,~{i-n-^ai+---)eP'^i~il-n-^ai+---)x^-2n^^3i[{djd,ai+---)^,xe''>)] 

^n{x):=-^ j^J^ isCC") (l + n-irfi+n-2d2 + ...)detP 

- Ix^)\)d£.dedV 
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Put h — ^, then we rewrite the above formula Dxih) := D„(a;). Dxih) is 
analytic with respect to h for fixed x. We Taylor expand Dxih) at /i = to obtain, 

D„(x) = D^{x) + O (-{x{l + x4)e-^')' 

7.4. Riemann surface. On a Riemann surface of area 1, we have 

(79) D„(x) = 4^xe-"' / e-l«l' |2|e|' 'X^\dC + O(-), 
or equivalently 

4 /"^ 1 

(80) D„(a;) = -xe"^ / e"'' |2r - x^l dr + O(-) 

TT Jo ' ' n 

8. Spherical ensemble: Proof of Theorem [2] 

In this section, we relate the expected density of critical points of an normal- 
ized random s„ G SH^{M,L'^) with the spherical Haar measure to the expected 
density D„ = DjJ" in the normalized Gaussian ensemble. We begin by relating 
Gaussian and spherical averages of density- valued random variables in Lemma [5] 
This relation is valid both for the critical point distribution and the value distribu- 
tion of ^ 

8.1. Relation between spherical and normalized Gaussian averages of 
density-valued random variables. 

Lemma 5. The spherical density and the 7^ densities of critical values are 
related by, 

^ni^) = w^2d„CD^(p-^a;)e-"V"-^c;p 
Proof. By ( 26 ) and Lemma [l] we have 

by definition of the normalization of the Haar spherical form on SH^{M, i")). We 
then change variables p — r^ . 

□ 

When a = c?„ , it follows from Lemma [5] that 

(81) T>i-{x) = i^d^"x2rf„-ijooj)^(^-i)g-d„.^p^rf„-i^^_ 

Changing variables to y = a;^, (81 ) is equivalent to 

(82) if„y-'*"+5d-'i"Df,"(2/i) = /:(p'^"-idg)(d„y), 
where C denotes the Laplace transform and we put 

(83) X„ = 27r'*"L.-; , dS(p) = Dn(p"')- 
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If we change variables — >■ p in (81 ) we also get, 

(84) /^„j/-'^" + 5d;^D^"(y5) = L{p''--iD^-id^s){y). 
In the last line, Drf{x) = f{r^^x) denotes the dilation operator. 

8.2. Spherical density for SU{2) polynomials. Before proving Theorem [2] it is 
helpful to do the calculations first for the case of SU{rn + 1) polynomials, where 
one has explicit formulae for ^pj^g formulae are simplest when m = 1 

and so we start with this case. We then follow the outline for the general case. 
For SU{2) polynomials, where d„ = n + 1, (|66]) and (|8T])-(|84]) imply that 

(85) i^„d„^2/-" (^2y-4— +8exp{-— y}j ^ Cip^^-^ D^-^d^s){y) 
We recall that 



(86) p-'e-^P = £lr„,, 



Hence the left side of (85) is Kndn ^ times 



^ (^21[i.oo](p)[ T{n) ] + «l[^^oo](p) J 



It follows that T)n'^^'''^\^/(^lP ^) — D,-i,„ equals 



We then put u = ^/d^p 2 to get 



Df'^^(')(ii) = if„d-iw (^21;o,^](^.)[u^d-i <V(l\"r - 4;r^ ^' rtn)"" ] + 81[o,6„] H ^^^P) ■ 
where fe„ = ^^g^^z^^- Recalling that (i„ = ?i + 1, we conclude that 
Df'^^(^)(.) = (21[„,^(u)[.2(n+ ~4^(i^^ . 

We further recall from (26) that W2(n+i) = '^ r{n+i) ^^'^ combine (n + l)^^T{n — 
l)-i ~ r{n + l)-i (resp. {n + l)-^T{n)-^ ~ r(n + 1)-^) to reach the final: 

Corollary 2. The spherical density for SU{2) polynomials of degree n is given 
exactly by 

Df'^^(^) = u {2u%o,V^ium - ^)'-2 - 4;^l[o,^(«)[(l - ^)"-i] + 81[o.b„](«)(l - ^7^)-^) 
It follows that 

lim„^oo Df'^^('V) = « (2"' - 4 + 8e-5"') e~"\ 
proving Theorem [2] in the SU{2) case. 
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8.3. Spherical density for SU{m + 1) polynomials. We now extend the proof 
to SU{m + 1) polynomials for general m. We do not simplify the integral for /,„ in 
that case but still may use the explicit n-dependence to confirm the main result. 

Let /3m — dimcS'(C'"), and let S'i(C™) denote the unit sphere in the space 
of complex symmetric matrices with respect to the usual inner product TtA*A. 
By Proposition [T] and since the determinant is homogeneous of order 2/3^ = 
dimR5(C™), we have 
(87) 



CmU^e 



7. 



det 



J,^ e-P |det [^\^^? - vp- 
^-fyhy^P-e-y 4 e-Py |det (^| VPcp - p-^) 



^drdu! 

'^f^^^-^dpduj 



^dpdiu 



■^yiy'^-e-y£{p'^-~'F^^-){y), 



where F„ 



\P) 



Isi 



det 



du! 



We further set = 3/3^ and recall (83). It follows from (84 1 that (with d^(p) 
Df'^^("+i)(p-^))) 



^ 2 ' 



■y 



^+V-e-yCip'^-~'F^'"'){y) = /:(p'^"-ii5 ,-.dg)(2/) 



We now simplify the left side using further identities for the Laplace transform. 
Denote the translation operator by Taf{t) = fit— a) and the Heaviside step function 
by H{t) = t^,. Also, denote the i/-fold primitive (fractional integral) by 



[x-ty 



-f{t)dt. 



We have (see e.g. [W] Theorem 8.1), 

£(i/(i-a)/(i-a))=e-"^£/(s) 



(89) 



We use the identities to simplify the left side of (88): 
(90) 

y-d„+a,„ + l£(^a,„-l^CP'")(^)g~y ^ y-'i"+'^'" + l/:Ti (if (p)p'''" - 1 ^^CP- ) (y) 

^ /:(/^.-i--..ri(ff(p)p'^"-iFCr'")(y). 



Combining ( 88 ) and ( 90 1 and uniqueness of the Laplace transform gives 

FT){P)=P'--Daz^d-M 



(91) 



n 2 



m 7-d„-2- 



"™d„Vi(i/(p)p" 
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We have, 
Therefore, 



1 

or equivalently with u ~ d?ip~2 , p — u^'^dn, 



(rf„-l-a,„) 



S,SU{m+l) 



Substituting the definition of we observe that 

(92) 

det(^^|VP^| 



= e-"V„°°e-*«^r"-i (4|det(|VPc.| 



- r 



dcj I (it. 



Also, 



hence 



r(d„ - 1 - am)(i° 



(d„ - 1) • • • (d„ - 1 - a„) 



r(d„) ~ r(d„). 



T{dn - 1 - a™) 



1 



Reversing the steps in (87) and comparing with (77), it follows that as n — > oo, 
(93) 



Cm „,2a™+l„-u" 



) dp 



~ Doo{u). 

8.4. Proof of Theorem [2} We now go through the general case, closely following 
the calculations in the special case of CP™. 
By Theorem [1] 

DjJ"(y5) = Doo(y^) + 0(^^), with Doo(2;') = /,„(y^)y^e-^. 



In determining the limit of it follows from (84) and (81) and the remainder 
estimate that we may drop the remainder term and then the calculation becomes 
almost identical to that of SU{m + 1) polynomials, with F!^^ (p) replaced by 



dio. 
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As in (87), we have 



(94) UWy) - ^y''-+'^e-yCip^--'F^)iy). 



By ([84]) ([94]) (with dg(p) = T>^{p-2))) we have 

(95) K^^dky'+-"^-''-{Cp---'Fm,){y)e-y = C{p'i"~i D^-icCi,)){y). 
Hence 

(96) ^-^K^I''-'^-'^-rrH{p)p'^-F„,{p) = p''--'^ D 

Repeating the calculation in the SU{m) case and setting u — V^nP^^ gives 
(97) 

We note that 

l[o.d„«--i](l - (^ + ^^)"' )^---°.--V"-if„(t) ^ exp(-(l +Ow')t""-^F™(i) 
montonically as n — )■ oo and so 

lim / " \l-^^^-*^)^"-''--V'"-iF„,(t)dt = e-"' /"exp(-pu2)^a,„-i^^ 







The rest of the calculation is the same as for CP™ and we get 

(98) lim 

9. Value distribution 

As a check on the results for the critical value distribution, we give the analogous 
calculation in this section of the much simpler expected value distribution, which 
is a probability measure on By the value distribution of (the modulus of) a 
section s„ € H'^{M,L") we mean the probability measure /Zs^ on defined on a 
test function ijj G C(]R+) by 

(99) Ps„ = i\Sn\h'^)*dV, {lP,{\Sn\h'^)*dV)= f ^bi\.Sn\h-)dV. 

JM 

As is well-known, the distribution is minus the derivative of the volume function, 

(100) y,,.(A) Vol{z e M : |s„(z)U. > A}. 

Fix z £ M and define the random variable pf^{x) = \ fniz)\e~"''^^^^^^ — |s„(z)|/in. 
For a test function tp e C(M+) we have (for any probability measure on sections) 

/"OO 

E„V'(pJ;) = / i^{x)f:{x)dx 
Jo 

where fn{x) is the pointise density of the distribution of p^. 

By the expected value distribution in the Gaussian ensemble we mean the mea- 
sures 

(101) A*n:=E„/i3„ 
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which is defined in the sense of distribution, 



(102) (V',E„(|S„U,.)*C«1^> - / ^n^KK{z)\h-)dV = En^{Pn)dV. 

JM JM 

To determine the expected value distribution it thus suffices determine E„V( I Sn(^) I 
or equivalently E„l/s^^(A). We note that /i„ is always an absolutly continuous mea- 
sure on K+ since s„ is analytic. Hence we may write it as fn{x)dx. 

We also define the expected value distribution in the spherical ensemble , 

(103) /xf = EfM,„, s„ e ^if"(A/,L"). 

It is also a continuous measure and we write /if = f^(x)dx. 

The random variable is the modulus of the Gaussian random variable 

which is only defined with a choice of local coordinates. But if V' is a radial test 
function on C, then 

E„V(pJ;) = E„V(p^), 
and the right side is globally well defined. We denote the expected density of by 
fniy), which is a density defined on C. We write it as f^{x,9) in polar coordinate 
y = xe*^. 

Lemma 6. For any probability measure on sections, the expected density of is 
given by the formula, 

(104) f^{x)= r f^{x,e)xde 

Ja 

9.1. Kac-Rice formula for the Gaussian value density. First we compare 
Gaussian value densities as a changes. We denote /" is the expected density under 
the Gaussian ensemble Q . The result is analogous to that of Lemma [2} 

Lemma 7. /"(x) = ai f^{aix). 

Hence it suffices to fix a ~ 1, and in the next Proposition we determine the 
Gaussian value distribution when a = 1 in Q ([5|. If we define the expected 
density of values , 

(105) El (^^^V(|snU")rf^^) =^ ^(a;)/i(a;)dx. 
Then we can prove, 

Proposition 2. The expected value density in the a = 1 ensemble is given by the 
formula, 

(106) flA^) = - I -^xe-^-^^^^y'^"dV{z), xe^+ 

TT Jm n„(z, Z) 

Proof. It follows from the Kac-Rice formula that 

(107) /^(y) = ^ exp {- {y, n„(z, z)-^y)) , 
since the covariance matrix of the random variable is 
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To obtain the density of values, we apply Lemma |6] and the polar coordinate y = 
xe*^ to (1071, then integrate the pointwise densities over M. 

□ 

The density is especially simple for SU{m+ 1) polynomials where n„(z, z) is the 
constant where Vm is the Fubini-Study volume of CP™. 

Corollary 3. For SU{m + 1) polynomials with a = 1 resp. dn and with volume 
Vm, the expected density of critical values is given by 

(108) /l^CP'"(x) = ^:.e-^-''"^-^ /^•«^^'" = 2Kna:e-^'"-^ 

As mentioned in the introduction, we usually set Vm — 1 for simplicity of nota- 
tion. 

9.2. Relation between spherical and Gaussian densities. The Gaussian value 
distribution has the same problem as the critical point distribution, namely the 
weighted repetition of sections. If we multiply s G 7?°(M, L") by r > then the 
volume function changes by Vrsi'^) — Vs{r~^X). Consequently, 

(109) frs^r-'fs{r-'s). 

The discussion is a word-for-word repetition of §8.1| Indeed, in the calculations 



we only used the relation ( 109 1 and the fact that the random variables take values 
in the densities on M. We therefore omit the proofs, but for clarity do state the 
analogous Lemmas. 

Exactly as in Lemma [5] we have 

Lemma 8. The spherical density and the 7^ densities of values are related by, 

f^{x) = if-ia'^" /;f (p-5x)e-"^p-5/'-idp. 

We combine Lemmas [7] and [H] to obtain 

Corollary 4. 

9.3. Spherical density of values: Proof of Theorem [3} In this section we 
prove Theorem [sj First we prove the statement for SU{m + 1) polynomials. 

Lemma 9. Let f^'^^ (u) be the density of values of SU{m + l) polynomials in the 
spherical ensemble with volume Vm ~ 1- Then 



hm /;f''^''"(u) = 2ue-"^ 

n— f 00 



Proof. The spherical density of values is determined from the Corollaries [3] and |4] 
(with Vm = 1) by changing variable y = ax'^, 

2K„^y-'-+'e-''-'y = £{p'--'^ F^){y). 

where we put 
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By (86) with a ~ d.^^ and = d„ — 1, we obtain 



Thus, 

With u = we get 



(110) 

The Lemma foUows 



/n - " r(d„-l) d„ -'-[d-\oo]l" "« i 



□ 



9.4. General Kahler manifolds. In the general case, we also have: 

Lemma 10. Let /^(u) be the density of values of random s G SH^{M,L"') in the 
spherical ensemble. Then 

lim (u) = 2ue""'. 



Proof. We combine Proposition[2]and Corollary|4] The exact formula of the Propo- 
sition shows that the general spherical density is asymptotic to the SU (m +1) case 
in dimenson m with a remainder term of order n^^(l + |x|'*)e~''"' . The remainder 



estimate is similar to but simpler than that for the critical point density in \ 7.3 
and is omitted. It follows that 

The rest of the calculation is then identical to the case of SU{m + 1) polynomials. 

□ 
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